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ABSTRACT 

X-ray and UV line emission in X-ray binaries can be accounted for by a hot corona. Such a corona 
forms through irradiation of the outer disk by radiation produced in the inner accretion flow. The same 
irradiation can produce a strong outflow from the disk at sufficiently large radii. Outflowing gas has 
been recently detected in several X-ray binaries via blue-shifted absorption lines. However, the causal 
connection between winds produced by irradiation and the blue-shifted absorption lines is problematic, 
particularly in the case of GRO J1655-40. Observations of this source imply wind densities about two 
orders of magnitude higher than theoretically predicted. This discrepancy does not mean that these ‘ther¬ 
mal disk-winds’ cannot explain blue-shifted absorption in other systems, nor that they are unimportant 
as a sink of matter. Motivated by the inevitability of thermal disk-winds and wealth of data taken with 
current observatories such as Chandra, XMM-Newton and Suzaku, as well as the future AstroH mission, 
we decided to investigate the requirements to produce very dense winds. Using physical arguments, hy- 
drodynamical simulations and absorption line calculations, we found that modification of the heating and 
cooling rates by a factor of a few results in an increase of the wind density of up to an order of magnitude 
and the wind velocity by a factor of about two. Therefore, the mass loss rate from the disk can be one, if 
not even two orders of magnitude higher than the accretion rate onto the central object. Such a high mass 
loss rate is expected to destabilize the disk and perhaps provides a mechanism for state change. 

Subject headings: accretion, accretion disks - hydrodynamics - methods: numerical - X-rays: binaries 


1. Introduction 

X-ray binaries (XRBs) are systems in which a 
compact object such as a black hole (BH) or neu¬ 
tron star (NS) accretes material from a secondary 
star. This material is believed to form an accretion 
disk surrounding the compact object, producing in¬ 
tense X-ray radiation via both blackbody radiation 
from the inner regions of disk itself, and Compton up- 
scattering of lower energy photons from a hot corona. 
These systems are further subdivided into two classes, 
based upon their properties. High-mass X-ray bina¬ 
ries have O or B class secondaries, are fairly steady 


X-ray sources and are thought to have evolved in-situ 
from a binary system. In contrast, low-mass X-ray 
binaries (LMXBs) have very faint secondaries, show 
dramatic changes in X-ray luminosity (X-ray bursts) 
and their evolutionary path is not clear. They may 
have evolved, through mass transfer and loss, from 
a situatio n where the donor star was m ore massive 
( Podsiadlowski, Rappaport & Pfahl 112002 ). or alterna¬ 
tively in dense regions like globular clusters it is possi¬ 
ble for a lone NS or BH to capture a low mass compan¬ 
ion (e.g. Verbunt & Lewin 20061) . LMXBs are highly 
variable systems, sometimes exhibiting an accretion 
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disk-like spectrum (referred to as the ‘high/soft’ state) 
and sometimes a power law spectrum (the ‘low/hard’ 
state). Sources in the low/hard state usually also show 
a radio jet which disa ppears as the source transitions 
into the high/soft state ( Fender. Belloni & Gallo 12004 , 


and refs therein) 

Accretion disks are predicted to have a flared 
geometry, meaning that the centrally generated X- 
rays will illuminate the surface of the disk, heating 


Shakura & Sunvaev 19731 

Begelman, McKee & Shields 

1983 

Ko & Kallman 1991 1 Raymond 1993; Woods et. al. 

1996 

Jimenez-Garate, Raymond & Liedahl 2002t 

Proga & Kallman 2002J) 

forming a hot disk atmo- 
al signature of this hot at- 
nd UV emission lines (e.g. 

sphere. The observation 
mosphere is both X-ray a 

Blair et. al. 1 1984; Raymond 1993; Jimenez-Garate et. al. 


2005 ) as the X-rays are reprocessed by not only the 
surface of t he disk, but also th e surface of the sec¬ 
ondary (e.g. Vrtilek et. al. 1990t) . This heating of the 
secondary has been proposed as a possible mecha¬ 
nism by which is mass transferred from the star to the 
accretion disk (e.g London. McCray & Auer 1 1981 


London & Flannery 1982). 

As observations of XRBs have improved, blue 
shifted absorption lines have been observed in more 
than_a dozen high/soft state NS and BH LMXBs 
(Diaz Trigo & Boirin 2013, and refs therein). These 
provide good evidence of the existence of outflow¬ 
ing material, most likely associated with the accre¬ 
tion disk, although the driving mechanism of these 
‘disk winds’ is a subject of ongoing discussion. A 
recent review of obser vatio ns of XRBs is given by 
Diaz Trigo_& Boirin (2013) and we summarize their 
data here. They cite outflow velocities between 400 
and 3000 km s _1 in 30% of NS systems, and 100 and 
1300 kms” 1 in 85% of BH systems. Photoioniza¬ 
tion modeling is usually employed to obtain estimates 
of the physical state of the absorbing gas, and such 
analysis gives a column density of between 4 x 10 22 
and 20 x 10 22 cm -2 and an ionization parameter of 
2.5 < log(£) < 4.5 for NS LMXBs. The ionization 
parameter is a measure of the ionization state of the 
gas, and we use the common definition 


L x 

nr- 


( 1 ) 


where Lx is the ionizing luminosity, n is the num¬ 
ber density of the gas, and r is the distance between 
the source of ionizing flux and the gas. BH LMXBs 


have a wider range of properties, with a column den¬ 
sity between 0.5 x 10 20 and 6 x 10 23 cuT 2 and an ion¬ 
ization parameter of 1.8 < log(£) < 6. The ion¬ 
ization parameter is degenerate in density and dis¬ 
tance, so in order to obtain information about the 
size/distance of the absorbers, it is necessary to break 
the degeneracy by measuring the density of the absorb¬ 
ing gas. T his has been done for the microquasar G RO 
.11655-4 0 dMiller et al lEoOhi iMiller et aklEoOsl but 
also see Netzer 2006) and these measurements sug¬ 


gest a relatively high density (n ^ 10 14 cm -1 ) which 
in turn implies a small radius. Those systems exhibit¬ 
ing absorptio n a ppear to be generally observed edge 


on (Ponti et. al. 12012) . which implies an equatorial 
geometry for the absorbing gas. As we will discuss 
later, this does not necessarily mean that any outflow is 
also equatorial - it can equally well be a bipolar flow, 
that exhibits stratification in physical properties such 
as density, ionization parameter or both. 

Possible mechanisms to drive disk winds are mag¬ 
netocentrifugal acceleration of g as g uided by mag- 
netic fields threading the disk (e.g. Blandford & Pavnel 


1982l:lEmmering. Blandford & Shlosman 1 1992 : Miller et al. 

2006al) . ra diation press u re acting on electrons or 
lines (e.g, Icke 1980t Shlosman & Vitello 19931: 


Shlosman, Vitello & Shaviv 19851 : Proga & Kallman 


2002) or thermal expansion of the hot disk atmosphere 
as a ‘thermal wind’ when the gas thermal vel ocity ex¬ 
ceeds the local escape velo city (IBegelman. McKee & Shields 
1983 : Woods et. all 1996 ). Whatever the mechanism, 
these winds are of great interest since they firstly pro¬ 
vide a way in which the XRB can interact with its 
surroundings, and secondly, if the mass flow is large 
enough, they coul d be the reason for t he observed 


state change (e.g. iNeilsen & Lee l2009t iPonti et. al 


2012 : King et. al. 20131) . The fact that these absorp 


tion features appear in high/so ft state sou rces but not 
in low/hard state sources (Diaz Trigo & Boirin 2013) 
is further evidence that they are linked to state change. 

In this work, we build upo n earlier simulations 
of thermal disk winds ( Luketic et al. 2010L hereafter 
L10). L10 modeled the launching of a wind in a sys¬ 
tem based on GRO J1655-40. In that work, a wind 
was launched but the velocity was too slow to account 
for the observed blue shifts of absorption lines and 
the density in the fastest parts of the wind was lower 
than the observed values for GRO J1655-40. Although 
these results suggest that thermal winds are unlikely 
to be the source of the absorbing material in that sys¬ 
tem, thermal driving is still an important mechanism 
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and deserving of further investigation. This is because 
it is almost certain to operate at some level in LMXBs, 
and even if it is not the principal source of the X-ray 
absorbing gas, it may well be important in the overall 
evolution of the system. For example, L10 showed that 
significant mass would be lost from the disk by thermal 
winds, and this mass loss is of the same order of mag¬ 
nitude as that which would be expected to d estab ilize 


the disk and perhaps drive state change (IShields et al. 
19861) . In addition, we cannot be sure if the wind in 


GRO J1655-40is typical or extreme, and new observa- 
tions which are likely come from the AstroH satellite 
( Done et. akll2014l Miller et aPl2014 , and refs therein) 
may provide more examples. We intend to investigate 
whether, with modifications to the heating and cooling 
rates, we can produce a thermally driven wind with 
physical properties more in line with current observa¬ 
tions, and provide a framework which may be of use 
in understanding future observations. 


The simulations we present here are all computed 
in the same way as described in L10, with three mod¬ 
ifications. Firstly, we consider only the optically thin 
case, so the radiation flux at any point in the simula¬ 
tion can be simply computed assuming a 1 /r 2 drop off 
from a centrally located source of X-rays. Secondly, 
we reduce the computational domain size from 20 Rjc 
in the original simulations to just 2 Rjc- Ric is the 
Compton radius, defined as the location in an accre¬ 
tion disk where the local isothermal sound speed at the 
Compton temperature, Tq, of the illuminating spectral 
energy distribution (SED) exceeds the local escape ve¬ 
locity. The Compton radius is therefore given by 


Ric - 


GM BH iun H 
k bTc 


( 2 ) 


where Mbh is the mass of the central object, equal to 
7M 0 for these simulations, // is the mean molecular 
mass which we set to 0.6, and other symbols have the 
usual meaning. The Compton temperature for the illu¬ 
minating SED in these simulations is Tq - 1.4x 10 7 K 
which gives a Compton radius of 4.8 x 10 11 cm. In all 
these simulations, as in L10, the disk is assumed to be 
flat and thin - defined via a density boundary condition 
at the midplane. 


The change in domain size is motivated by prelim¬ 
inary investigations which show that the acceleration 
zone for the wind was located inside 0.1 Ric, in line 


with 

irevious work (e.g. 

Begelman, McKee & Shields 

1983 

Woods et. al. 1996; Proga & Kallman 2002 

). In 


essence, this is because the most efficient acceleration 


of gas via thermal expansion occurs as the gas is heated 
past the lower equilibrium temperature on the thermal 
equilibrium curve (see Figure[2]and related discussion 
in the next section). This occurs at fairly low values 
of f, which is where the gas is densest, i.e. close to 
the central object. The gas then enters an unstable 
heating zone where the next stable temperature is over 
an order of magnitude higher. Rapid heating occurs 
resulting in rapid acceleration as the gas expands. This 
change to the domain size means that we better simu¬ 
late the densest parts of the wind where absorption is 
most likely to occur. We neglect the outer parts of the 
disk which means that our calculated mass loss rates 
are lower limits. 

Finally, but central to this project, we will mod¬ 
ify the heating and cooling rates assumed in the heat¬ 
ing term in the thermodynamic equations. We exam¬ 
ine several cases, each with different rates, and each 
representing a modification to the thermal equilibrium 
curve. This changes the way in which the gas passes 
through the ‘heating’ region, and we will see that this 
can have a profound impact on the velocity, density 
and hence mass loss rate of the wind. 

In the next section, we briefly discuss our method¬ 
ology. We then discuss the details of the flows pro¬ 
duced by the different heating/cooling parameters. Fi¬ 
nally, we discuss the relevance of our results to the on¬ 
going discussion of thermal wind in XRBs. 

2. Method 

Our simulations are based upon the simulations pre¬ 
sented in L10 (run C8). Fig. Q] shows a rerun of that 
simulation, but computed upon a smaller grid, concen¬ 
trating on the inner parts of the flow where the acceler¬ 
ation occurs. The wind in this simulation, and all oth¬ 
ers presented here is shown at a time of 220 000s which 
represents 25 sound crossing times ( Ric/c s ), a time 
that is sufficient for the wind to have settled down to a 
steady state. Note that this timescale, about 2.5 days, 
also gives us an estimate of the lag between a change 
in the luminosity or SED of the central source of ion¬ 
izing radiation and the response of the wind through a 
change in structure. Of course changes in the ioniza¬ 
tion state of the wind would likely occur more quickly. 

There are several important features of this model 
that can be seen in the figure. Firstly, we see a dense, 
turbulent ‘corona’ at to < 0.3R/c, z < 0.2Rjc where 
gas flow is severely inhibited by gravity and reaches 
steady state only in a time averaged sense. The stream- 
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lines starting just outside this region can escape, how¬ 
ever it is still time dependent, hence the streamlines 
show ‘kinks’. 


Secondly, outside the inner corona, the flow starts 
off moving vertically, as gas expands away from the 
accretion disk which lies along the lj axis. The stream¬ 
lines bend outwards, partly due to the pressure gradient 
and partly as a result of conservation of angular mo¬ 
mentum which means that the centrifugal force acting 
on the rotating gas is not balanced by gravity. They are 
self similar in this region. A fast, fairly dense flow is 
produced which escapes at angles less than about 45°. 
Finally, we see a fast, low density infall at polar angles. 


Although the L10 simulation did produce a fast 
fairly dense wind, as mentioned in the introduction the 
density was 2-3 orders of magnitude lower than that 
inferred from observations of GRO J1655-40. Even 
so, the total mass loss rate through the wind was seven 
times the accretion rate. We use the same version of 


the hy drodynamics code ZEUS-2DlSto ne & Norman 
d 1992 ) extended by Proga, Stone & Kallmanl ( 200(1 ) to 
carry out 2.5D simulations of the flow. In this code, ra¬ 
diative heating and cooling of the gas is computed us¬ 


ing the parametrized rates first described by [Blondin 


19940 . This parameterization includes Compton heat¬ 
ing and cooling, photoionization heating, recombina¬ 
tion cooling, bremsstrahlung cooling and collisional 
(line) cooling as functions of temperature T and ion¬ 
ization parameter A The ionization parameter is de¬ 
fined as in Equation Q] The number density of the 
gas is related to the density of the gas by n = 
where the mean molecular weight p is set to 0.6. Opti¬ 
cal depth effects are not considered in the simulations 
here, so the factor of Lx/r 2 is related to the flux, re¬ 
duced only by distance effects. 


Compton heating and cooling is given by 


G Compton — 8.9 x 10- 36 A7 x - AT) (ergs s' 1 cm 3 ) (3) 


where Tx is the temperature of the illuminating power 
law spectrum (set to 5.6xl0 7 K). Photoionization heat¬ 
ing and recombinational cooling are subsumed into 
one term, referred to as “X-ray heating/cooling”, given 
by 

G x = 1.5 x 10~ 21 £ 1/4 r~ 1/2 (l - T/T x ) (ergs s" 1 cm 3 ) 

(4) 

whilst bremsstrahlung cooling is parametrized by 

L b = 3.3 x 10“ 27 r 1/2 (ergs s _1 cm 3 ). (5) 


Finally, line cooling is given by 

L, = [l.7 x 1CT 18 exp (~T L /T)C l T~ l/2 + 10 24 ] 5 

( 6 ) 

where 7/ has the units of temperature and parametrizes 
the line cooling. It is set to 1.3 xlO 5 K. The 5 parameter 
allows one to reduce the effectiveness of line cooling 
due to opacity effects. In an optically thin plasma, 6 is 
set to 1. The units of L/ are the same as for the other 
rates. 

We are interested in investigating whether simple 
changes to the heating and cooling rates, thereby mod¬ 
ifying the thermal equilibrium curve, can increase the 
velocity and density of the wind to better match obser¬ 
vations. To modify these rates, we apply pre-factors 
to each of the mechanisms and so the equation for the 
net cooling rate £ (ergs s - g 1 ), which appears in the 
energy conservation equation, becomes 

P-C - n ~(A c Gcompton + A X G X - A/L/ - A b L b ). (7) 


The first six lines of Table [j] gives the values of 
pre-factors, Tx and Lx for the original L10 simu¬ 
lation (run C8, denoted A) and six further simula¬ 
tions. The line cooling pre-factor effectively replaces 
the 5 parameter and therefore represents a measure 
of line opacity. Modifications of Ax, the photoion¬ 
ization/recombination rate pre-factor can be justified 
by a change to the illuminating SED or metallicity of 
the gas. Calculations of the precise nature of the con¬ 
nection between these parameters and Ax are beyond 
the scope of this work, our values are not intended to 
represent a particular case, rather we adjust them to 
produce the desired thermal equilibrium curve. We 
change A b somewhat arbitrarily to make gas thermally 
stable everywhere. Our aim here was to investigate the 
relationship between thermal instability (TI) and effi¬ 
cient acceleration. The upper and lower stable temper¬ 
atures remain the same and so we isolate the effect of 
instability with this experiment. 


The quantity £ C oid,max is the value of the ionization 
param eter when the flow becomes thermally unstable. 
Field I ( 19651) demonstrated that in a non-dynamic al 


flow, subjected to isobaric perturbations, thermody¬ 
namic instability results when [6£/6T] p > 0. This 
is also where the gradient d\n(T)/d\n(^) becomes 
greater than 1. Table|T|also gives T eq , the equilibrium 
temperature expected for £ C oid,max- It can be shown that 
if line cooling is balanced by X-ray heating on the cool 
branch of the stability curve, then this temperature is 
expected to be 4/5 T b . This is 104 000K, and look- 
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Fig. 1.— The density structure of model A (see Table |T|). Overplotted are streamlines (grey lines), the 80° and 
60° sightlines (dashed lines) and arrows showing the velocity field. Also shown is the location of the M=1 contour 
(black line). 
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ing at the values in Table |T] we see that the actual 
values are very close to this. We also include values 
of E C oid,max, which is the ratio of radiation pressure to 
gas pressure when the flow becomes thermally unsta¬ 
ble^ This is_given by H = Fj on /nki,Tc = f/ Ank^Tc 
( Krolik, McKee & Tartar! 1981), where the tempera¬ 
ture T is set to the equilibrium temperature at the onset 
of instability. 


To produce comparable simulations, we use a den¬ 
sity boundary condition to ensure that the ionization 
parameter is equal to £, co id,max at the Compton radius. 
The density at the midplane at the start of our simula¬ 
tions is defined by the equation 


pO) = po 



where po is given by 


Po 


L x m H fi 


£coldjnaxR 


2 ' 

1C 


( 8 ) 


(9) 


This is given in the table along with Ric for each case. 
Since the density in the disk is proportional to 1/r 2 , 
this means that f is a constant in the disk plane, and 
it is the same for all runs except Ah in which it is ten 
times bigger. 

The hydrodynamic calculations are carried out in a 
spherical polar coordinate system, running from 0 to 
90° in angle, and from R,„ to R out in the radial direc¬ 
tion. The zone spacing increases with radius, such that 
drk+i/drk = 1.05 giving finer discretization in the in¬ 
ner parts of the flow. The zone spacing reduces with 
increasing angle dOk+i/dQk = 0.95 giving more resolu¬ 
tion close to the disk. These parameters, together with 
the number of points used in the two dimensions are 
given in Table [Q 

The solid lines on Figure [2] show the thermal equi¬ 
librium temperature predicted for the different cases 
plotted against £JT, which is equal to the ratio be¬ 
tween the radiation and gas pressures. In an outflow, 
the gas pressure cannot increase along a streamline. 
This means that £ IT must always increase]], and so as 
a parcel of gas is heated, starting at the bottom left of 
the graph, it will follow these curves until the gradi¬ 
ent becomes negative (the onset of TI, if present). This 


'Note that this is in fact only true in the case of weak or 
zero magnetic fields. If magnetic fields are strong, then the 
gas pressure can increase and the gas can move to the left on 
the equilibriu m curve (E mmering, Bland ford & Shlosm an 119921 : 

iBottorff. Korista & ShlosmanT^OOol) 



Fig. 2.— Thermal equilibrium curves for the seven 
cases considered. The solid curves show the equi¬ 
librium temperature (T) vs £/7\ For each case, the 
crosses represent actual data from the 60° sightline 
and the circles are from the 80° sightline. The size of 
the symbol shows the radial distance along the sight¬ 
line, with larger symbols representing the largest dis¬ 
tances. The triangles show the point at which the heat¬ 
ing curve becomes unstable. 


location on the curve represents the maximum temper¬ 
ature of the cold branch. The triangle symbol on the 
graphs shows this point. At this point, one expects the 
gas to quickly heat up to the upper equilibrium tem¬ 
perature. This maximizes the rate of energy transfer 
between the radiation field and the gas, driving expan¬ 
sion and hence acceleration. Thus it is in the unstable 
zone where the most ‘efficient’ acceleration of the gas 
takes place in order to form an outflow. 

The behavior of the different cases is best explained 
in the context of the shape of these thermal equilibrium 
curves. Case A and Ah have identical thermal equilib¬ 
rium curves, becoming unstable at the same value of £ 
and reaching the same upper equilibrium temperature 
(set by the balance of Compton heating and cooling). 
The difference in luminosity between the two cases 
has no effect on the shape of the equilibrium curves, 
however as we will see in the next section, it does 
affect how the gas is heated. Given the same unsta¬ 
ble zone, but enhanced radiation density, one would 
e xpect case Ah to produce_a faster, denser flow (see 
Begelman, McKee & Shields"!! 19831) . 

Case B has reduced line cooling and enhanced 
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Prefactors 

A 

Ah 

B 

C 

D 

E 

F 

A, 

1.0 

1.0 

0.2 

1.0 

1.0 

1.0 

0.076 

Ac 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

A b 

1.0 

1.0 

1.0 

3.9 

1.0 

1.0 

1.0 

A x 

1.0 

1.0 

4.0 

1.0 

1.0 

1.0 

1.0 

Physical Parameters 

T x ( 10 6 K) 

56 

56 

56 

56 

0.80 

230 

295 

L x (33 x 10 37 ergs s _1 ) 

1 

10 

1 

1 

1 

1 

1 

log ( £ cold,max ) 

2.10 

2.10 

0.91 

N/A 

N/A 

2.07 

1.2 

lo§(. ^ cold,max ) 

1.33 

1.33 

0.17 

N/A 

N/A 

1.32 

0.43 

T eq (tcold,max)( 10 3 K) 

111 

111 

106 

N/A 

N/A 

109 

113 

po(10' 12 gcm“ 3 ) 

1.14 

11.4 

17.4 

1.14 

1.14 

22.4 

281 

/?/c(10 10 cm) 

48.2 

48.2 

48.2 

48.2 

3380 

11.5 

9.15 

Grid parameters 

cm) 

2.4 

2.4 

2.4 

2.4 

2.4 

1.2 

1.2 

^m fl .v(10 10 cm) 

96 

96 

96 

96 

96 

96 

96 

R ratio 

1.05 

1.05 

1.05 

1.05 

1.05 

1.05 

1.05 

Nr 

80 

80 

80 

80 

80 

80 

100 

@min 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

@max 

90.0 

90.0 

90.0 

90.0 

90.0 

90.0 

90.0 

® ratio 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

N e 

100 

100 

100 

100 

100 

100 

100 

Wind properties 

V r (max blueshif ted )/100 km s _1 

4.47 

6.76 

6.78 

4.28 

2.11 

14.5 

16.3 

V r (p > lel2, max blueshif ted )/ 100 km s _I 

1.18 

1.98 

3.26 

0.203 

0.186 

5.03 

2.51 

rin (60° sightline)(x 10 22 cm -2 ) 

5.71 

53.0 

74.2 

1.01 

0.281 

13.3 

15.0 

riH (80° sightline)(x 10 22 cm -2 ) 

46.3 

476 

441 

22.6 

24.1 

16.6 

255 

M wind,disk ( M acc ) 

3.72 

3.61 

41.5 

0.878 

0.723 

4.20 

25.7 

M w ind,outer(Macc) 

3.45 

3.36 

38.3 

0.638 

0.758 

3.85 

24.9 


Table 1: The heating and cooling parameters adopted in the simulations, and some key parameters of the resulting 
winds. 
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photoionization heating. This means that the gas 
becomes unstable at a lower Changes to the 
SED have been shown to have this effect (e.g., see 
Bottorff, Korista & Shlosman 200(1 who use detailed 


photoionization calculations), so our approach is not 
unreasonable. The distance to the radiation source will 
be largely unchanged, so low £ equates to higher den¬ 
sity and so one would expect that the resulting outflow 
would be denser. 


Case C and D are both attempts to remove the un¬ 
stable zone entirely. Case C has strongly enhanced 
bremsstrahlung cooling, through an increase in A*. We 
did this, somewhat arbitrarily, so that there is no TI for 
10 5 < T < 10 6 but the equilibrium temperatures at the 
low and high end are the same as in A, Ab and B. Our 
goal here was to isolate the importance of TI to accel¬ 
eration of the gas. Case D achieves the same thing, by 
reducing the X-ray temperature. This reduces both the 
initial photoionization heating rate, and also the upper 
equilibrium temperature. 

Finally, cases E and F represent solutions with two 
unstable zones. In case E, we simply increase the X- 
ray temperature. This does indeed produce two unsta¬ 
ble zones, however the second unstable zone is at a 
lower pressure than the first. In case F, we increase the 
heating rate on the lower branch in a similar way to 
case B, by reducing line cooling. This shifts the lower 
unstable zone to lower £ but leaves the upper unstable 
zone unchanged. The aim for these two cases is to see 
if one can get a faster wind by extending the accelera¬ 
tion zone. It is also interesting to see if gas ‘collects’ 
at the stable zone between the two unstable ones. 


3. Simulation results 

We present the results of our simulations in the con¬ 
text of the thermal equilibrium curves. This allows 
us to see how the hydrodynamics of the winds affects 
the thermal balance, and thus give insight into how the 
winds are accelerating. The symbols, plotted over the 
solid lines, on Figure [2] show the relationship between 
the temperature and ionization parameter divided by 
that temperature in a range of cells along two sight¬ 
lines. We can therefore see if the equilibrium temper¬ 
ature is reached in the simulations. 

We note that there are no points on the lower sta¬ 
ble branch of the stability curves for any of the cases. 
This is by design, since points along the lower branch 
are essentially in the disk, and merely exhibit turbu¬ 
lent motions. We have carried out detailed resolution 


studies using ID simulations which demonstrated that 
the behavior of the wind in the regions sampled by our 
sightlines is not dependent on resolving the transition 
from cool stable branch to unstable zone. Our limited 
2D resolution study confirmed this. 

Turning our attention to individual cases, we first 
examine cases A and Ah. The increase in luminosity 
has had the desired effect, in that all of the points for 
case Ah lie on the upper branch of the unstable zone, 
whilst those of case A lie below the curve. Adiabatic 
cooling as the hot gas expands holds the temperature 
below the expected stable temperature. This indicates 
that more energy has been transferred to the gas in the 
higher luminosity case. Note however that the points 
from the two simulations occupy remarkably similar 
locations in T - £ space, given that Ah has an order 
of magnitude higher luminosity. This is because the 
density of the wind has increased accordingly, giving a 
very similar value of the ionization parameter for both 
flows. 

Case B has, as expected, produced hotter gas at 
lower However, as with case A, the gas is cooler 
than the upper stable temperature. This suggests that 
with a higher luminosity, one could obtain even hot¬ 
ter gas. Interestingly, the 80° points are all very close, 
indicating that we are sampling points at the same rel¬ 
ative distance along each streamline that the sightline 
intercepts. 

Cases C and D, where there is no unstable zone, the 
points tend to lie along the thermal equilibrium curves 
with some exceptions. The innermost 60°sightline 
data of case C are cooler than expected if radiative pro¬ 
cesses dominated the heating/cooling. This is because 
adiabatic expansion is acting as an additional cooling 
mechanism. This simulation fails to launch a wind in 
the current simulation domain. Case D also fails to 
launch a strong wind and the polar infall which is al¬ 
ways seen in these simulations extends to much larger 
inclination angles. Compression of the gas by this flow 
at low radii produces some gas that is heated above the 
upper stable temperature. This occurs for gas with an 
ionization parameter greater than about 10, and is not 
shown on Figure [2] 

Case E shows two new effects. Firstly, despite 
there being two unstable zones, the gas is unable 
to access the upper unstable zone since this lies at 
higher pressure. Therefore the gas jumps straight 
from the lower unstable point towards the upper stable 
branch. Adiabatic cooling prevents the gas from reach¬ 
ing the upper branch and the existence of gas in a for- 
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mally prohibited part of the thermal equilibrium graph 
demonstrates that hydrodynamic effects are an impor¬ 
tant consideration in the thermal balance of this type of 
wind. It is often assum ed th at gas w ill avoid unstable 
zones (e.g. Chakravorv. Lee & NeilserJ 20131 and refs 
therein), providing an explanation of why some ion¬ 
ization states are not seen in observations. This sim¬ 
ulation demonstrates that the situation may be more 
complex. 

Finally, case F does produce points close to the sta¬ 
ble zone between the two instabilities since that part of 
the curve is now physically accessible to the flow. 

In the lower section of Table [T] we provide some 
of the physical properties of the simulated outflows. 
First, we list the maximum radial velocity seen in each 
of the simulations. It is of the same order of mag¬ 
nitude (a few hundreds of kilometers per second) as 
the blue shifted absorption lines seen in LM XBs (e.g. 


Diaz Trigo & Boirin 120131 : iNeilsen 2013, and refer¬ 


ences therein), however the velocity of gas is not the 
only important factor. 


The density is also important since only dense gas 
will produce observable absorption. The maximum ve¬ 
locity in regions with a particle density greater than 
10 12 cm -1 is much smaller. This is simply because we 
are probing regions deeper into the wind, where the 
flow is still accelerating, and we see that the two sta¬ 
ble cases do not produce fast enough gas at high den¬ 
sity. Cases B, E and F do show fast moving gas with 
density above the threshold density. In cases E and F, 
the fast, dense gas is limited to a very narrow range of 
angles, in case E it is within 3 degrees of the disk and 
would therefore probably not be observable. In case E, 
the fastest dense gas is at small radii around 6 -15°. 
Whilst this could in principle be observed, the small 
angular range over which an absorption feature would 
appear would likely mean it would be transient. 


Another important physical parameter that can be 
derived from observations is the total column density. 
This is between about IQ 2 0 - 10 23 cm) 2 for all kinds 
of LMXBs ( Diaz Trigo & Boirin 20131) . Our simula¬ 
tions produce column densities of the right order of 
magnitude for equatorial sight lines, and indeed the 
80°sightline would be Compton thick in case Ah and 
B. 


Finally, we give the mass loss rate, both leaving 
the disk ( M W ind,disk ) and leaving the computational do¬ 
main (Mwind,outer)- These rates are calculated directly 
from the simulation results, using M = £ Af ,v where 


A is the area represented by a cell, either on the disk 
or at the outer boundary, p is the density of the cell, 
and v is either the vertical velocity for disk cells, or 
the radial velocity for cells at the outer boundary. 
The summation is carried out over all relevant cells. 
We report these values in terms of the accretion rate, 
M acc = 4.4 x 10 17 g s -1 (assuming an efficiency of 
8.3%). In L10, their version of model A produced an 
outflow of about 7 times the accretion rate whereas we 
see an outflow rate of about half that much. This is 
simply because our simulation has a much reduced ra¬ 
dial extent compared to the L10 run. Increasing the 
luminosity (model Ah) increases the mass loss rate, 
but the ratio of mass loss to accretion rate remains the 
same. Cases B and F produce significantly higher mass 
l oss r ates, in excess of the threshold of 15M flec that 


Shields e t al. (119860 demonstrated could induce oscil¬ 


lations in an accretion disk. It should be noted that the 
winds are emerging from the disk far outside the radius 
where most of the radiation is produced. Almost all of 
the ionizing radiation from a thin disk in a system like 
this one is produced within 100 gravitational radii of 
the centre. By comparison, the induced Compton ra¬ 
dius is about half a million gravitational radii. There¬ 
fore, at least for the purposes of these simulations, the 
model of a point-like, unvarying source of radiation at 
the centre of the simulation is not necessarily invali¬ 
dated by the prediction of large mass losses from the 
outer parts of the disk, at least over the timescale of the 
simulation. 


4. Synthetic absorption line profiles 

It is useful to produce synthetic line profiles for our 
simulations, in order to get some idea of whether the 
outflows could, in principle, produce the absorption 
features observed in XRBs. Computation of the ion¬ 
ization state and level populations of the gas is beyond 
the scope of this work, and we calculate the absorp¬ 
tion using the simplified scheme discussed below. This 
scheme takes account of thermal line broadening, and 
the doppler shifting due to the bulk flow of the wind, 
however we do not model line emission here. 

The opacity due to a resonance line (uncorrected for 
stimulated emission) is given by 

a{v) = ~r~n j Z?i 2 0(v) (10) 

47T 

We arbitrarily set Bn = 1 and use the hydrogen num- 
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ber density for n i. Therefore, our line opacity becomes 
hv 

a v =—n H <p(v). (11) 

An 


For the line shape 0(v) we use a gaussian line profile 
of the form 


< P(y) = 



2nkgT 


exp 


me 2 (v - vo) 2 'l 
2 k B Tv 2 ) 


( 12 ) 


For each radial cell i, we compute the line opacity as 
a function of frequency, and then doppler shift that 
line profile to take account of the bulk flow velocity. 
We then obtain the total, frequency dependent optical 
depth by summing up the opacity at each frequency 
due to each cell of radial thickness dr. 


i=outer 

t (v) = ^ an(v)dr. (13) 

i=inner 

This sum is computed from the innermost radial cell 
to the outermost, thus making the implicit assumption 
that the continuum source is point-like, and located at 
the origin. A measure of the absorption profile of a 
generic line is then computed 

F(v) = e~ T(v \ (14) 


Since no attempt is made to compute the density of any 
given ionic species, these spectra are in no way accu¬ 
rate representations of what we would expect to ob¬ 
serve for a given system. Rather, they are just a means 
of comparing runs, since each spectrum is calculated 
in a consistent way. 

The upper panels of Figure [3] show the results of 
the line profile calculation for the 60° (left panel) and 
80° (right panel) sightlines. We immediately notice 
that the 80° features are much deeper than those seen 
at 60°. This is simply due to the higher density at the 
base of the wind. This means that the ionization pa¬ 
rameter is also generally lower in the 80°sightline, and 
so it is likely that different species would be observed 
at the two angles. 

A general feature of most of the spectra is an ab¬ 
sorption feature close to zero velocity, and a second 
feature at a blueshift that varies from model to model. 
Although absorption features are seen at zero velocity 
in observations, it is fair to say that the very strong fea¬ 
tures we predict here are not commonly observed. We 
know however that by including line emission these 
absorption features would be weakened. We are more 


interested in the second absorption feature which ap¬ 
pears at velocities between about -150 km s _ and 
-800 km s^ 1 . 

Looking in detail at each of the cases now, we first 
see that the increased luminosity of case Ah over case 
A was partially successful in that the absorption is 
stronger. This is of course due to the increased den¬ 
sity of the outflow. This can be seen clearly in the 
lower left panel, where the cells in case A have density 
less than 10 12 cm -3 (represented by the cross symbol) 
whereas the same cells in case Ah have density greater 
than 10 12 cm 3 (circles). However, increase in den¬ 
sity was not our only aim, we also wanted to increase 
the velocity at which the absorption was seen. In this 
we have been only partially successful, with the blue 
shifted absorption feature shifting to only very slightly 
higher velocities. 

Case B is perhaps the most successful new model 
with blue shifted absorption features seen in both the 
60°and 80°sightlines. The blue shifted absorption fea¬ 
ture at 80° is deep and its width is only due to thermal 
broadening. This is because all of the cells produc¬ 
ing the feature are at nearly the same velocity. This 
is in turn because photons flying along this sightline 
encounter gas in a very similar physical regime at all 
radii showing that, close to the disk at least, the gas 
is flowing along highly self similar streamlines in this 
case. Since the ionization parameter is set to be the 
same at all radii at the mid plane, all gas will start in 
the same physical state. In case B, the physical state of 
the gas has evolved similarly at all radii, ‘remember¬ 
ing’ its initial conditions, up to the 80°point. In con¬ 
trast, by the time it has moved up to the 60°sightline, 
that ‘memory’ of the starting state has been lost, and 
gas at different radii is in different physical conditions. 
Thus the absorption is produced by cells at a range of 
velocities and so the feature is much shallower but very 
broad. 

As already discussed above, case C and D fail to 
produce fast outflows. This is clearly shown in the 
lower panels, where all the cells from these simula¬ 
tions are clustered around zero velocity. These two 
cases produce relatively narrow (the temperature of the 
gas is lower) features at zero velocity. 

As shown in Table [j] cases E and F do both pro¬ 
duce dense, fast moving material. However, the mate¬ 
rial in case E is very close to the accretion disk, and 
is missed by both sight lines shown here. Appearing 
at angles greater than 87°, it is unlikely that it would 
be observable in any case. Case F does produce fast 
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60 ° 80 ° 



Fig. 3.— Simulated spectra for the 60°(left hand upper panel) and 80°(right hand upper panel) sight lines. Lower 
panels show ionization parameter vs radial velocity for all cells in the two sight lines. The size of the symbol represents 
the distance from the central mass ( small symbols are for small distances). Crosses represent cells with a number 
density of less than 10 12 cm~ 3 and filled circles show cells with a density above this threshold. Negative velocities 
represent motion away from the centre of the simulation, and thus blue shifted absorption. Note that in the bottom 
right hand panel, the black and blue points overlie one another. 


11 

















material at lower 6 and this is seen in the spectra as 
absorption around -200 km s _1 . Faster gas does exist 
in the simulation, but only over a very narrow range of 
angles. 


5. Discussion 


Our aim was to see if simple, physically motivated 
changes to the heating and cooling balance of the 
thermal wind simulated in L10 could produce a wind 
model that was more in line with observations. There 
are three main observational measures, the line veloc¬ 
ities, the column density and the density of the line 
forming region. 


The wind model described by L10 failed to pro¬ 
duce any gas with velocities greater than 100 km s _ 
with a density greater than 10 12 cm~ 3 and we come to 
a similar conclusion. Since L10 were trying to repli¬ 
cate the properties of the outflow observed in GRO 
J1655-40, which seems to have a ve ry high density 
of 5 x 10 15 cm -3 ( Miller et al. 2006a ) the conclusion 
was that the model failed in that aim. Follow up 
work reduced the d e nsity estimate to aroun d 10 l4 cm 3 
( Miller et al. 2008; Kallman et. al. 20091) but this is 
still well above the densities seen in the L10 model. 
This failure of a thermal wind model to replicate the 
velocity and density seen in the observations sug¬ 
gests that for GRO J1655-40 at least, a thermal wind 
seems an unlikely source of the observed absorbing 
gas. Nonetheless, the predicted wind produces a col¬ 
umn density in line with observations, and a mass loss 
rate in excess of the accretion rate. 


The enhanced luminosity version of the L10 model, 
case Ah, does increase the velocity and perhaps more 
importantly the density of the wind. A radial veloc¬ 
ity of 200 km s 1 of gas with a number density greater 
than 10 12 cuT 3 is still slow and less dense compared to 
the measurements of GRO J1655-40 mentioned above, 
but it is not unreasonable to think that this wind would 
produce observable features, and further work to char¬ 
acterize the ionization state of this wind allowing cal¬ 
culation of detailed spectra would be worthwhile. The 
ionization parameter for the fastest moving parts of the 
wind has a narrow range, centered on log £ ~ 4, and is 
certainly similar to that inferred for the absorbing gas 
in many systems. 

Case B provides the best illustration of how sim¬ 
ple changes to the heating and cooling rates in a sim¬ 
ulation can affect the velocity and density of the re¬ 
sulting wind. We reduced line cooling by a factor of 


5, and increased the photoionization heating rate by a 
factor of 4. Both of these changes can be broadly jus¬ 
tified, by the effects of line optical depth in the first 
case, and changes to SED and gas metallicity in the 
second. This simple change made the gas thermally 
unstable at £ one order of magnitude lower. The radial 
location of the unstable gas is largely unchanged, so 
the change in H means that denser gas is accelerated, 
producing a denser and faster wind. Although the den¬ 
sity and velocity is only a little higher than case Ah, 
much more interesting is the huge increase in mass loss 
rate, now 40 times the accretion rate (even though we 
only simulate a relativel y small domain). Th is is al¬ 
most 3 times the rate that lShields et alJ (119861) showed 
would induce instabilities in the disk. Therefore, even 
if thermal winds are unable to reproduce the observed 
line absorption seen in XRBs, they may well provide a 
mechanism for XRB state change and so searching for 
an observational signature is a worthwhile exercise. 

Another interesting result from these simulations is 
that there is a gas with ‘forbidden’ values of <(', i.e. 
from the second, hotter, unstable zone of the stabil- 
ity cu rve in c ases E a nd F. It has been suggested (e.g. 


Chakravorv. Lee & Neilsen ; 2013) that species which 


are expected to have peak abundances in gas with such 
forbidden ionization parameters would not be seen in 
observations. Whilst our results do not necessarily dis¬ 
prove such assumptions, they do illustrate that hydro- 
dynamic effects (i.e. adiabatic expansion) make the 
situation more complex. 

It is often assumed that disk winds in XRBs are 
equatorial, because absorption is_ seen preferentially 
in sources which exhibit dips Ponti et. all ( 2012 ). We 
find that the wind is in fact bipolar, but the outflow 
is highly stratified with the high density region of 
the wind near the disk. Therefore, in our simula¬ 
tion, absorption is only seen for equatorial sight lines 
even though the wind flows out over a relatively wide 
range of angles. Similar results were seen in L10 and 


Giustini & Proga (120121) who computed line profiles 


based on LlO’s as well as other disk wind simulations. 
This stratification is also important with respect to 
the observed variability of X-ray absorption in XRBs 
( Neilsen 20131 and refs therein). This is very well 
illustrated by cases E and F, both of which produce ab¬ 
sorbing gas in narrow angular ranges. If the illuminat¬ 
ing spectrum is variable, the angle at which particular 
species would be seen could change - thereby making 
the absorption lines associated with those species van¬ 
ish. The wind could remain strong in this case, but 
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would need to be detected in different species. 

Secondly, when we compare cases A and Ah, which 
differ only in luminosity, we see that the density of 
the wind solution changes significantly giving rise to 
a very similar ionization parameter. This also has rel¬ 
evance to studies of variable sources, where increases 
in luminosity are sometimes called upon to explain in¬ 
creases in £ and hence the disappe arance of some fea¬ 


tures (e.g. Diaz Trigo et al. 20141) . Our results show 


that it may be overly simplistic to assume the density 
remains constant in such cases, and a more detailed 
investigation is required, taking into account how the 
wind responds to the increase in luminosity. 

6. Future work 


The simulations we have presented here use a 
simplified heating and cooling scheme, which per¬ 
mits swift exploration of parameter space. In ad¬ 
dition, radiative transfer through the wind is treated 
in the optically thin limit. Previous detailed analy- 


sis of s_uch hydrody namic models (Sim et al. 201(A 


Higginbottom et. al.l 2014) have shown that a more 
thorough treatment of radiative transfer including scat¬ 
tering can have a significant effect. We therefore plan 
to run such simulations on the more successful models 
from this work (i.e. those that produced high velocity, 
dense flows). Not only will this work give more in¬ 
formation regarding the validity of our modified heat¬ 
ing/cooling rates, but it will also produce detailed ion¬ 
ization data for the wind and spectra. It will also pro¬ 
vide information regarding line emission contribution. 
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